Adiabatic-antiadiabatic crossover in a spin-Peierls chain 
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We consider an XXZ spin-1/2 chain coupled to optical phonons with non-zero frequency ljq. In 
the adiabatic limit (small ujo), the chain is expected to spontaneously dimerize and open a spin 
gap, while the phonons become static. In the antiadiabatic limit (large u)o), phonons are expected 
to give rise to frustration, so that dimerization and formation of spin-gap are obtained only when 
the spin-phonon interaction is large enough. We study this crossover using bosonization technique. 
The effective action is solved both by the Self Consistent Harmonic Approximation (SCHA) and by 
Renormalization Group (RG) approach starting from a bosonized description. The SCHA allows 
to analyze the low frequency regime and determine the coupling constant associated with the spin- 
Peierls transition. However, it fails to describe the SU(2) invariant limit. This limit is tackled by 
the RG. Three regimes are found. For ujq <C Aj,, where As is the gap in the static limit wo ~* 0, the 
system is in the adiabatic regime, and the gap remains of order As. For ujo > As, the system enters 
the antiadiabatic regime, and the gap decreases rapidly as ujo increases. Finally, for ujq > ujbkt, 
where ujbkt is an increasing function of the spin phonon coupling, the spin gap vanishes via a 
Berezinskii-Kosterlitz-Thouless transition. Our results are discussed in relation with numerical and 
experimental studies of spin-Peierls systems. 



I. INTRODUCTION 

The properties of the spin-Peierls (SP) state in 
quasi-one-dimensional materials has attracted consider- 
able attention over the last decades since its discov- 
ery in the organic compounds of the TTF and TCNQ 
seriesiSi^, and more recently in the inorganic compound 
CuGeOg^. In analogy to the Peierls instability in quasi- 
one-dimensional metals^, a spin-chain undergoes a SP 
transition by dimcrizing into an alternating pattern of 
weak and strong bond o^i^i^ with the magnetic energy gain 
compensating the energy loss from the lattice deforma- 
tion. 

Although this physical picture gives a good qualita- 
tive understanding of the SP phenomenon, the real SP 
transition is in fact much more complicated to describe. 
In particular, the above picture of SP transition is only 
valid in the adiabatic regime in which the frequency of 
the phonons is negligible compared to the magnetic en- 
ergy scales in the system, such as the spin gap or the ex- 
change interaction J. The validity of this approximation 
is clearly dependent on the system at hand. Recently, it 
was pointed out that the difference between CuGeOs and 
the other SP compounds consists in the high energy of 
the optical phonons involved in the transition, which is of 
the order of the exchange integral J"'^^'^^'^^). Another fea- 
ture that distinguishes CuGeOa from the other organic 
SP compounds is that no softening of the phonon modes 



is observed near the transition. All these findings stem 
for the fact that an adiabatic treatment of the phonon 
subsystem^ is inadequate to describe the SP transition 
in CuGeOa and an appropriate treatment of phonons in 
the antiadiabatic regime is requiredi^ii^. 

Unfortunately, not many analytical methods to study 
the system of coupled spin and phonons in the full fre- 
quency range are available. The main difficulty relies 
in the fact that when the phonon frequency becomes 
comparable to the energy gap in the spin-excitation 
spectrum , one is entering a quantum regime in which 
quantum fluctuations completely impregnate the ground 
state. This is why many of the known studies in- 
volving dynamical phonons rely on numerical methods 



such as exact diagonalization (ED) 
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strong cou- 



pling expansions.i-, density matrix renormalization group 
(DMRG)22*2iiS2i2^ or Quantum Monte Carlo (QMC) 
simulationsS^iSSiSSiSLS^. From the analytical point of 
view various approaches have been developed, but they 
work well either in adiabatic or in the antiadiabatic 
regime. In the former case, most approaches are based 
on the mean-field approximation2i22i2£. In the latter case 
various perturbation studies were performed to derive an 
effective spin Hamiltonian^iiSS. Another approach was 
developed, based on integrating out the phonon modes^^ 
in order to map the model onto the Gross-Neveu model24 
for which various exact results are available^^^^^. Other 
approaches are based on the fiow-equation methodic that 
works well in the antiadiabatic regime, or on the unitary 
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transformation method for the XY spin chain^i. 

Since various compounds are rather close to the bor- 
der between the two regime^^, it would be thus highly 
desirable to have a good method to tackle the adiabatic- 
antiadiabatic crossover. In this paper we provide such 
a method. We combine the renormalization group (RG) 
method and the self-consistent-harmonic-approximation 
(SCHA) to study the adiabatic-antiadiabatic crossover 
in a spin-1/2 Heisenberg chain coupled to dynamical 
phonons. A previous attempt to use the RG to study this 
adiabatic-antiadiabatic crossover has been published^^. 
We will comment on the differences of the two ap- 
proaches. 

The plan of the paper is as follows: in Section UTI we 
introduce the model of spin chain coupled to dynami- 
cal phonons, and write it in the continuum limit using 
the bosonization representation. In Section IIIII we de- 
scribe a variational approach, inspired from the self con- 
sistent harmonic approximation (SCHA), and use it to 
describe the crossover. In the adiabatic regime we find 
an expression for the spin-Peierls gap consistent with the 
mean- field treatment of Cross and Fisher^. In the an- 
tiadiabatic regime the gap is essentially described by a 
sine-Gordon model. In section Hvl we study the crossover 
using a renormalization group (RG) method. The RG 
is specially well adapted around the Heisenberg isotropic 
point. This allows to extract the phase diagram as a 
function of the strength of the clcctron-phonon coupling 
and the phonon frequency. In section [3 we discuss the 
findings of the two methods, both relative to each other 
and in connection with experiments. Conclusions can be 
found in Section lVll and some technical details have been 
put in the appendices. 



II. MODEL AND CONTINUUM LIMIT 

As a simple model which describes a SP system in the 
following we consider an antiferromagnetic spin-1/2 chain 
coupled to a set of Einstein oscillators. 



with 



H — Hs + Hp + Hgp, 



Vn , m 2„2 

2 m 



(1) 

(2) 
(3) 



where Sn are spin 1/2 operators, J > 0, [qn,Pn'] = i5n,n' ■ 
The quantity muj"^ — is the stiffness of the Einstein 
phonon. The interaction of spins with phonons can be 
modeled by: 



(4) 



The coupling to optical phonons described by 10} is ad- 
equate for the CuGeOa since it would correspond to a 



side group effect by Germanium atoms as discussed in 
Refs. . 39^40 . Acoustic phonons could of course be treated 
in a very similar way, but one would have to replace Qn 
with {qn+i — Qn)- Note that some authors^^-^® prefer to 
diagonalize the phonon Hamiltonian in ^ using the bo- 
son operator b — ^ mjjjQ/2q + ip/ ^2mujQ, and write the 
interaction g'Y^^{b\^ + bn)Sn ■ S„+i. It is obvious that 
one has g = gj \J2mL0Q. This remark will be useful when 
we will compare the results of the different approaches. 
The adiabatic limit is — *■ 0, m — > oo with fcg fixed. 
In that limit, the phonons become classical, i.e. the q„'s 
commute with the Hamiltonian (Q, and one can simply 
minimize the ground state energy with respect to their 
expectation value. In that limit, the results of Ref. ^i^are 
recovered. The opposite antiadiabatic limit is — > oo. 
In that limit, one can integrate out the phonons, and 
one is left with the Hamiltonian of a frustrated spin 1/2 



chai: 



For a frustration large enough: 



41.42.43 



for 



large enough spin-phonon coupling, a spontaneous dimer- 
ization of the spins takes place, and the system presents 
a spin gap. Our purpose is to provide a unified treatment 
of these two limits. 

To solve the spin-phonon problem, we use first the 
well-known Jordan- Wigner transformation to express the 
spin-operators in terms of spinless fermions. Thus the 
Hamiltonian iJ^ becomes: 

Hf = -iI]k+iCn + h.c.] + 

n 

V^(c];+ic„+i - i)(4c„ - i), (5) 

n 

with t ~ J /2 and V = J . The spin-phonon Hamiltonian 
Hsp is transformed into: 



^(ci+ic„ -l-4c„+i) -t- 

(4+ic„+i - ^)(4c„ - i) 



(6) 



We now proceed in the standard way to take the contin- 
uum limit (see e.g. Ref. i44 chap. 6). In the continuous 
approximation, © generates a coupling between the lat- 
tice deformation (phonon mode) and the q = 2kp = tt 
component of the charge density, p(2kp,x). 

In order to get a continuous description we separate 
fast and slow components of the phonon field and sim- 
ilarly for the fermion fields, so we get the interaction^ 



Hfp = i / dx[q(x)p{2kp,x) ~ h.c] 



(7) 



We now use the boson representation of one dimensional 
fermion operators. In this representation the Hamilto- 
nian Hf becomes 

Hf = ^j dxuK{^n{x)f + ^iyd,{x)f (8) 
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where the field (f>{x) is related to the density of fermiong^ 
and n(x')] = i6{x — x'). We have u — ^Ja, with 

a the lattice spacing, K = 1/2, g„ = 17(2: — na) and 
we have kept only the most relevant terms. Changing 
the parameter K allows one to explore the more general 
case of XXZ spin chains with an easy plane anisotropy^. 
The long wavelength part of the fcrmion density is 
Pqr^o{x) — — iV0(a;) whereas the higher Fourier com- 

P2fe,(x) = — (J)'cos(2(/.(x)), (9) 

where we have specialized to an isotropic spin chain. The 
prefactor in Eq. Q has been shown in Refi^ to yield a 
good agreement of the gap calculated within bosonization 
and numerical calculations''^. The Matsubara action for 
the phonon field has a standard quadratic form: 

Sp = ^ j dx dT[{drqf + Lolq\ (10) 

I 



where p = m/a\s the mass density of the optical phonon 
mode, and q(x = na) = qn is the lattice deformation 
field. In the first approximation (0) we have neglected 
the fact that the phonon disperses. It can be shown that 
the dispersion along the chain leads to insignificant cor- 
rections. On the contrary if the phonons are three dimen- 
sional, i.e. if uj disperses with the transverse momentum 
then significant changes can occur. Indeed in that case, 
since the phonon are three dimensional they couple the 
different spin chains and can induce a three dimensional 
transition at low temperatures. We will come back to the 
case of three dimensional phonons in Section Ivl 



Since the total action is quadratic in the phonon 
fields, we can integrate them out to obtain the follow- 
ing bosonized action with a retarded interaction between 
the electronic densities: 



S = dx 



dr 
2ttK 



r r0 rl3 

i dx i dr dr' cos24>{x,T)T>i^g_f3{T — t') cos2(j){x,T'), 
J Jo Jo 



r 



where "D^^^fjij) is the propagator of an Einstein phonon 
of frequency corresponding to the action (|10(l 



V^,,p{r-r')^{q{T)q{r')) 
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]r^r'\ , 2cOsh(tJo(T~ rQ) 



(12) 



The action l|ll|) fully describes a one dimensional spin 
chain coupled to phonons, and does not rely on adiabatic 
or antiadiabatic limit. However, one has to note that 
because of the cutoff, the action Hll|) is valid only for 
Wq «C u/ a. For higher value of the phonon frequency ujq, 
the phonon propagator H12|l must be replaced with 5{t ~ 
t'). In that case the action Hll() is simply the continuum 
action of a frustrated spin chain*^-, in agreement with the 
canonical transformation approachi^*^. 

The action Hll|l is of course impossible to solve exactly. 
In order to obtain the physical properties of the system 
we will analyze it using two different techniques in the 
next sections. The first technique is a self consistent ap- 
proximation. Such an approximation, will be very useful 
to define the various phases of the system as well as the 
relevant parameters. As any variational approximation, 
although it can be very efficient in describing the various 
phases it can only describe the transitions between these 
phases approximately. So in order to study the critical 
points we use a renormalization group method, building 
on the knowledge of the relevant parameters extracted 
from the SCHA. 

III. SELF-CONSISTENT HARMONIC 
APPROXIMATION 



ational method^°i^^i^^i^^. The idea is that the action l(TT)) 
would be classically minimized by <^ = 0. One can thus 
expect that the physics will be dominated by small devi- 
ations around this minimum and approximate the action 
((TT|l by a quadratic action. We thus consider as a trial 
action 

So = ■^^^G^'^{q,UJn)4>*{q,(^n)(l>{q,l^n) (13) 



We have to find the propagator G{q,ujn) so that (|13l) is 
the best approximation for Hll|l . For that we define the 
variational free energy 



^var. — Fq + {S — 5o)c 



(14) 



where: 



Fo 



i_lnZo = -i^lnG(c.„), (15) 



I3L 



D(j>e 



-som 



(16) 



and (. . .)o represents an average with respect to the ac- 
tion Sq. 

The second term in the action (|ll|l can be rewritten as 



To study the action Hll|l we apply first the Self Consis- 
tent Harmonic Approximation (SCHA) or gaussian vari- 
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2(7ra)2pajQ 



/3 /■/? 

dx / dr I dr'2?^o,^(r - T')[cos(2(/)(a;,r) + 2(/)(x,t')) +cos(20(x,r) - 2(/)(a;,r'))]. (17) 











r 



Given that the phonon propagator (|12|) decays for large 
time difTerence (r — r'), one can see from (|17ll that the 
cosine of the sum is roughly equivalent to ~ cos(40(a::, r)) 
and can be responsible for the opening of a gap in the 
spectrum, while the cosine of the difference is ~ {t — 
T')^(VT0(a;, t))^ and thus will modify the quadratic part 
of the action. 

In the following we consider the gapless (A = 0) and 
the gapful (A ^ 0) case separately, at zero temperature. 
The gapless case is interesting in connection with systems 
of electrons at an incommensurate filling interacting with 



phonons^^iS^i^. In these systems, the term cos[2(j){x, t) + 
2(j){x, t')] does not appear in H17() . leading to A = 0. The 
spin-Peierls problem corresponds to the half-filling case 
for the fermions. 



A. Incommensurate case 



Using (|13|l the variational free energy is given by: 



= ^{lGo\'l,^)-G-\q,u;)]G{q,u;)-\nG{q,u;)} 



2/ dr^e-"l-l(cos(20(O,r)- 0(0,0))), (18) 



where GQ^(q,w) — -i^iuq^ + ^) and we used (fT^ for 
(3 ~ oo. Introducing the propagator for the field 0, 
G(x,r) = (0(x, t)0(O, 0)), we can rewrite: 



(cos 20(0, r) cos 20(0,0)) = 



(19) 



Using this expression we obtain: 



£*fe-".Mexp{-4/|i/|^G(,,.,|l-c„.,„ 



^r)] . (20) 



Minimizing the action 120|l with respect to G(q,uj) and where S is the self-energy, we get: 
using the fact that G^^{q,Lj) = GQ^{q,iLj) — T,{q,iLj), 



SFyar „ 1 f dq f did 

^G(^=°^2y 2;^y 2^[^('^'") 



(7ra)2pw2 



rf^^^e-°'-l(l-cos(c.r))exp|-4 f ^G(g,a;)[l - cos(^t)]|]. (21) 
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As is obvious from the above equation, the low frequency 
behavior of G((7, w) is similar to the one of Gq^ {q,uj), and 
corresponds to a variational action of the form 



So 



dx 



dr 
2ttK 



(22) 



thus G ^{q,uj) = ■^{uq'^ + %)■ equation above 

we thus have: 



dq f du) [1 — cos(a-'T)] 



ttK 



K 



211 J 27r {uq^ + ^) 



Huct), (23) 



where = ^, is a frequency cutoff. Using (|23() we can 
write our variational equation as: 



T.{q,Lo) + 



9 



^ (u;or) ^_^„l^l (l-cos(^T)) 
2 (wct) 



2K 



0. 

(24) 

If we confine to an expansion up to order in (1 — 
cos(ijJt)), as requested by the analytical behavior of the 
Green's function for cj — > 0, we obtain: 



T.{q,uj) 



9 



(7ra)2pcj2 J 

^2 



(^^cT) 



2K 



9 



2K 



(25) 



where T is the gamma function. As we see the integral is 
convergent when K < 3/2. Going back to the definition 
of the self-energy we have: 



( 1 



1 



V 2-KuK 2'KuK 



(26) 



Equating (|25|l with (|26|l , and using the fact that 
we obtain the following value of the parameter K: 



if' 



2Kg^ ( u \ 
TTupojQ \ aujQ J 



2-2K 



r(3-2is:)]. (27) 



Expanding around K , we obtain the renormalized value 
of if: 



K'l- 



TTupujQ \ awo J 



(3-2K)]. (28) 



One thus recovers a Luttinger liquid but with a renor- 
malized value of the Luttinger parameter K. H28|) implies 
that K < K. A similar result can be obtained via the 
renormalization group analysis (see the next section). In 
a RG analysis the result would correspond to in- 
tegrating the RG equation for the coupling constant g 
assuming that K is not renormalized and then comput- 
ing the lowest order correction to K with the renormal- 
ized coupling constant. Our method thus reproduces in 
a cr ude way t he renormalization of K downwards. As in 
Ref. I55l56l57l we find that the tendency of the system to 
form charge density waves is increased. 

B. Commensurate case 



In the commensurate case the derivation of the varia- 
tional free energy from (|14I21|) remains the same in the 
commensurate case. Let's rewrite (|20ll in slightly differ- 
ent way: 



J^.a. = i/ ^/ ^[G^\q,u;)G{q,cj)-\nG{q,u;)] 



^{'KQfpUjl 2 



iexp|-4y"^ J ^G{q,Lj)[l + cos{LJT)] 



Minimizing this action with respect to G{q, lu) we obtain the following expression for the self-energy 
E(,, .) . £ ..^f e-M ((1 - cos(..)) exp {-4 / I / ^G(., .)[! - cos(..)] 



(29) 



+ (l + cos(tJT))exp|-4 J ^ J ^G{q,Lj)[l + cos{ujT)]^^ (30) 
I 

As is obvious from the above equation H3()() . E(g, cj) is in expansion for the self-energy: 

fact independent of q. Moreover, we can use the following i 

E(9,c.) = -_(A2+7c.2). (31) 

TTUK 
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In Eq. (|31|l . the variational parameter A stands for the 
gap caused by the commensurabihty, and the variational 
parameter 7 stands for the renormalization of the bare 
Luttinger exponent K. Such a restricted ansatz is justi- 
fied by the fact that higher powers of lo in Yj{uj) are asso- 
ciated with irrelevant operators in the action, whereas 
A and 7 correspond respectively to a relevant and a 
marginal operator. Keeping only A amounts to neglect 
any renormalization of K by the spin-phonon interaction. 

The self-energy leads to a Green's function 

G{q,w): 



where Kq is the Bessel function. The corresponding vari- 
ational action is 



5*0 = I dx 



dT 
2nK 



1 ?7 

u 4 



,(34) 



(32) 



where A is the mass term. The integral of the Green's 
function is: 



where u/^ = A is the gap and u/K 
{dx(j>)'^ is generated from lfT7|) . 



u/K as no term 



Equating the coefficient of H31() with that coming from 
the expansion for small oj of (|3()|l , we obtain the following 
two equations: 



29' 



_ ^ 

ttK {naYpuiQ 

1_ 2ff^ 
■KuK 



, 2 



dq 
2^ 



duj 
2^ 

G{q,uj) [1— cos(c(Jt)] 



G'(g, w)(l -I- coswr) 



(35) 

(36) 



Using that — the l.h.s. of H3()|) can also be rewritten 



as: 



1 



1 



TTUK 



2ttuK 2TruK 



(37) 



The two self-consistent equations (I35|l and (|36|l can be 
solved analytically in the antiadiabatic limit (luq 3> A). 
Using H33|l . and after a straightforward but lengthy cal- 
culation we obtain: 



Kg^ f u \ 
TTupuiQ \aujo J 



2-2K 



r(3-2if)]. (38) 



which is the same change of K than in H28|) . The system 
also develops a gap given by 



A = ^ 



nupujQ \auJo J 



(39) 



As we can see from (|38|l . for K > 1/2, we can have 
K < 1/2, so that H39(l can still lead to a gap provided 
that g is large enough. Combining the two equations l|38|l 
and H39|l . we finally have: 



,2-4K 



T{3-2K) 
T(TT2K)^ 



(40) 



The SGHA thus correctly describes the formation of a 
gap in the antiadiabatic limit. As for the incommen- 
surate case the SGHA captures part of the effects of the 
renormalization of the parameters. Note that the SGHA, 
as any variational method is efficient in capturing the na- 
ture of the ordered phases, but in order to determine the 
nature of transition one needs the full RG analysis. Such 
an analysis will be discussed in Sec. IIVI 



C. Adiabatic-antiadiabatic crossover in the SCHA 



Using the SGHA we are now in a position to describe 
the crossover from adiabatic to antiadiabatic regime. We 
will assume that we are far from the point K = 1/2 and in 
fact that we have K < 1/2. For the spin chain this would 
correspond in being in the Ising limit. In that regime, we 
can neglect the renormalization of K and take K = K in 
our variational action. The variational free energy H29() 
can be written: 



7 



F = Fn 



:G(0,0)- 



9 



e 'a 



r 



2K 



g-4G(0,r) _^ g4G(0,T) 



(41) 



where G'(0,t) is given by Eq. (j33l. Using the expansion 
for the Bessel function we obtain the following approxi- 
mate expression; 



G(0,r) = ifMT>2e"'^^, 



if UT < 2e-'<i, 



442) 



where 7 is the Euler-Mascheroni constant^^. Using the 
above expression for the Green's function, we are able to 
calculate the variational free energy H41|) , and minimizing 
it with respect to ^ we obtain the following variational 
equation: 



e 'a 



2K 



+ 



e 'a 



r 



AK 



\aLOQ 



2K 



7(1 



u 



(43) 



where 7(-, •) is the incomplete Gamma function^. 

Two interesting limits in equation l|43|l must be dis- 
cussed. If ^i^fii ^ 0, one is in the adiabatic limit, whereas 
the antiadiabatic one corresponds to — > 00. In the 
adiabatic limit, one sees that the term on the right hand 
side of (|43ll reduces to a contribution ^ (a/S,)'^^ , the term 
depending on the incomplete Gamma function being zero 
in that limit. As a result, the Cross-Fisher prediction for 
the gapS, 



A 



2k, 



(44) 



is recovered. In the antiadiabatic limit, the exponen- 
tial term in l|43|l disappears, and the incomplete Gamma 
function can be replaced by a Gamma function, leading 
to the result for the gap we have found in Sec. IIII Bl in 
()39|l . In this limit, the gap can be understood as result- 
ing from a cos 40 interaction induced by integrating out 
the phonon modes. 



To perform a general study for any loq^/u, we rewrite 
(gSIl for £_ as: 



(f) 



2-2K 



^ \2i^oi J 

In terms of the gap, this equation reads: 



2-2K 



where: 



^2-2K 



fix) 



AKg^ 

TTUpUjQ \ 2 



2K 



+ (xf)^^7(l + 2if,2^ 



2K 



(45) 



(46) 



(47) 



The graph of the function f{x) is represented on Fig. ^ antiadiabatic regime is easily observed, with the two lim- 
In this figure, the crossover from the adiabatic to the iting forms of the gap given respectively by Eq. H44|) and 
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FIG. 1: The graph of the function f{x) (solid line) defined in 
Eq. (133 with K = I. Two regimes are visible: For A 3> wo, 
f(x) ~ x^^^^ (dashed curve). In that regime, the gap is 
given by the adiabatic formula Eq. 1441 . For A <^ uo, f{x) ~ 
j.2-4,K ^(jQ^^g(j curve) and the gap is given by the antiadiabatic 
formula Eq. The crossover regime is observed for 0.3 < 

A/loo < 3. 



(|39|l . The SCHA allows to get the full interpolating func- 
tion between the two regimes, and thus to obtain pre- 
cisely the crossover scale. We obtain that the limit be- 
tween the adiabatic and the antiadiabatic regime is given 
hy ujQ A and not by loq ^ J. This point will be further 
discussed in the forthcoming section Hvl 

The SCHA also yields the expectation value of the 
nearest neighbor correlations S„ • S„+i, as it is propor- 
tional to (— l)"(cos 2(/>). One finds: 



{Sn ■ Sn+l) 



K 



For loq <^ A, i.e. in the adiabatic regime, one has: 



{Sn ■ Sn+l) 



In the antiadiabatic regime, for K < 1/2, wc find: 



(48) 



(49) 



{Sn • Sn+l) 



g' ( u \ 



2K 



(50) 



of the quadratic part by the phonon coupling term. Such 
a renormalization of the parameter K is of course spe- 
cially crucial to take into account precisely close the 
isotropic Heisenberg point K = 1/2. In this section, 
we thus apply an RG method to analyze the adiabatic- 
antiadiabatic crossover. 

Attempts to an RG analysis of such a problem or 
of directly related fermionic problems have been de- 
scribed in the literature. In particular an RG analysis 
was performedJ^ at T = based on a previous work 
on spinful fermions coupled to phononsSiffiiSi. In this 
work, the interaction of the spinful fermions with the 
electrons is viewed as a retarded backscattering interac- 
tion. However, although this description is appropriate 
for fermions, in the case of the spin chain it neglects the 
fact that the staggered dimer operator gives rise to more 
relevant interactions than current-current ones. As a re- 
sult, this fermionic description underestimates the size 
of the dimerization gap. Our analysis, directly based on 
the boson representation of the spin chain does not suffer 
from such a limitation. In addition to providing us with 
a better description of the dimerization gap, the use of 
the boson representation also allows us to tackle the case 
of a finite frequency lvq and nonzero temperature. 

Another closely related problem is the one of fermions 
in a random potential^, which has an action quite simi- 
lar to (|11|) but with a constant V. One could be tempted 
to simply reuse the RG equations derived for this sys- 
tem. However, here, the situation is more subtle. In lO 
for uJo/T , V{t) T. Therefore, we see that the 
rescaling of the temperature is going to modify the RG 
equations with respect to the case of disordered fermions. 
Moreover, for ujq/T oo, I?(r) = wo/2e~"ol^l. As a re- 
sult, the limit of T — > is delicate to handle properly. 
In particular, the definition of the spin-phonon coupling 
constant becomes ambiguous in this limit. 

However the variational analysis performed in the pre- 
vious section allows us to build the correct RG proce- 
dure. First the variational approach shows that in order 
to obtain the correct results it is important to first per- 
form the calculation of the ground state free energy for 
< T <C A, where A is the spin-Peierls gap, and then 
take the limit of T ^ 0. Second it gives us that the 
proper dimensionless coupling constant measuring the 
strength of the electron-phonon interaction is: 



IV. RG ANALYSIS 



TTUpUjQ 



(51) 



As we have discussed in the previous section, the 
SCHA describes only approximately the renormalization 



We now proceed with the RG. We start from the following 
action: 
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dx I drcosAcj). (52) 





The cos(40) operator is the marginal operator needed to 
describe an spin-isotropic spin chain. The derivation of 
the equations is given in Appendix O They read: 
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2' (53c) 



(53d) 



These RG equations are conveniently expressed using the 
coupling constant G defined in Eq. I|51|l . At this one loop 
order, we find no corrections to the phonon frequency as 
can be seen in (|53dp . However, we expect such correc- 
tions to be obtained in a higher loop order calculation. 



A. Anisotropic case 

Since the action also describes spinless fermions 
coupled to phonons, our equations have similarities with 
the RG equations that have been derived for the elec- 
tron phonon problem^^i^SiSiSSiS. There are however im- 
portant differences. First for a spin chain the equiva- 
lent fermionic band is automatically half-filled (in the 
absence of an external magnetic field). Thus in ad- 
dition to the standard terms that were considered for 
the electron-phonon problem with incommensurate fill- 
ing, one has here to take into account the marginal umk- 
lapp operator cos(40) as in Ref. Second, in the 

electron-phonon problem a different coupling constant is 
usedSS*^, namely ^ Q^^oi^^ guch definition ap- 
pears natural when looking at the RG equations H53a() 



and (|53b(l since Ygp seems to be the amount by which K 
is renormalized in the limit T = 0. However, such defi- 
nition would be at odds with the calculations performed 
with the SCHA. In fact, the integral dl^V^„^A^) = 1 
for all LOo- As a result, if wc neglect in (|53a|) . and the 
renormalization of K in l|53c|l . wc find the following ap- 
proximate RG equation for G and K: 



G{1) = G(0)e(2-2^)', 



(54) 



d 
Jl 



{K-') = G(0)e(2-2^)'^e'expf-i^e')>(55) 

u \ u / 



and by a variable change to 1^ = ^^^e', we easily obtain 
that 



K-'{oo)-K-\0) = G(0) r(3-2i^). (56) 



2-2K 



Wo a 



This equation is easily understood: ujo gives an energy 
cutoff that stops the RG fiow of K induced by G at an 
energy scale of order ojo = u/ae''' . We note that it is 
identical to the SCHA result (^5)) . We thus see that at 
that scale, K is renormalized by an amount proportional 
to G{1*) and not G{1)^^^^^ as a result of the exponential 
factor in 155|) . This confirms that the right coupling con- 
stant in this theory is G and not Gujo{l)oi/u. In Ref. [gJ 
the same prescription was used to define the coupling 
constant whereas in Ref. the incorrect rescaling of 
Ref. I55I wa s used. As a result, we expect the conclusions 
of Ref. I33 to be incorrect in the adiabatic regime. 

Till now, we have assumed that at the scale I* = 
ln(u/(Q;cL'o)), the coupling constant G{r*) ^ 1. If this 
assumption breaks down, since the coupling constant 
G{1) = e(2-2-ft^)'G(o) one finds a gap 



A = ^G(0)i/(2-2K) 



(57) 



This gap is in agreement with the SCHA result and with 
the mean field theory treatment of Cross and Fisher^ . For 
K < 1/2, in the antiadiabatic limit uo 3> u/a, we know 
from the SCHA that the phonons can generate a relevant 
perturbation cos 40 and thus induce a gap^i. This effect 
is also captured in the RG by (j53bp . This can be seen by 
a two step renormalization procedure. In the first step, 
for I < I* =ln — a term g 1 is induced by the RG flow. 
This term is found to be of order: 



TTU \ aujo ) 

7(2if+l,l)-7(2i^+l,— ) 

u 



(58) 



Since ujo ^ w/a, we can actually neglect ^(2K -I- 1, ^^) 
in Eq. jSS)). For I > I*, 'D^^i^i^{a/u) 0, and we can 
drop G from the RG equations. We then have a simple 
Kosterlitz Thouless RG flow which leads to a gap of the 
form: 



G(0) 



\aLOo 



2K 



7(2X+1,1) 



(59) 
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This gap is in agreement with the SCHA prediction in the 
antiadiabatic hmit (|38|l . Therefore, we see that SCHA 
and RG methods agree perfectly, far from the isotropic 
point, once the proper couphng constant is used in the 
RG. 

Using our RG equations we can now study the SU(2) 
invariant hmit for which the SCHA cannot be used, due 
to importance at that point of the marginaUy irrelevant 
operator cos(40). 



B. SU(2) invariant case 



In the isotropic limit, we have 



K 



2 V 2ttuJ 



(60) 



This ensures that, in the absence of spin-phonon cou- 
pling, the flow will renormalize to the fixed point K* — 
1/2 and g*l — 0. It is then easily seen that the equations 
(|53al53b|) reduce to a single equation for 2/ = f^- This 
leads to the following RG flow: 



dy 
dl 



y 



G{l)—ee 



dG / 3 , ^ 



(61) 
(62) 



These RG equations allow for the full interpolation be- 
tween the adiabatic and antiadiabatic limit. 

The simple analysis of the previous section showed that 
the gap should behave as A = ^G{Q) in the adiabatic 
limit. For the isotropic case, using (|61I62|) . we obtain 
logarithmic corrections to the gap A = —G\ InGj^'^/^ re- 
sulting from the marginal flow of y{l). These logarithmic 
corrections (for details see the Appendix ID|) are identical 



to those obtained by incorporating the logarithmic cor- 
rections to the gap of the dimerized spin 1/2 chain^SiiMl 
into the Cross-Fisher mean field theory. This confirms 
that G is the right coupling constant to study the forma- 
tion of the spin Peierls gap in the adiabatic limit. On the 
other hand, as discussed in the previous section, in the 
antiadiabatic limit, it is the fiow of y{l) that determines 
whether or not the gap is formed. To analyze the fiow in 
the antiadiabatic regime, we can use the approximation 
G{1) — G(0)e' i.e. we neglect the logarithmic corrections 
to the flow of G. We have checked that this approxima- 
tion leads to a good agreement with the numerical study 
of the RG flow using the fourth order Runge-Kutta algo- 
rithm. Using the previous approximation, the RG flow 
()61|I - H62|) can be reduced to a Ricatti differential equation 
(cf. Appendix IeJ leading to the following dependence of 
the gap on G: 



A 



■ Woe' 



exp 



uG{Q) 



(63) 



for the case of y(0) = 
a gap exists only if: 



0. When y{Q) < 0, it is found that 



uGjO) 
2uJoa 



> 



|y(o)l 



1 



l2/(0)|ln(^ 



(64) 



The physical content of this equation is transparent. At 
the scale I* such that wqc' = u/a, G{1*) is equal to the 
l.h.s. of the inequality whereas |?/(Z*)| is equal to the r.h.s. 
of the inequality. The gap can form only if the renormal- 
ized spin-phonon interaction is stronger than the renor- 
malized marginal coupling at the energy scale loq. This 
is in agreement with the two step RG approach^ of the 
preceding section. When the condition H64|) is satisfied, 
the gap behaves as: 



A = wne-(^-^) 



exp 



mG(0) 
2LUQa 



yjo) 



(65) 



l-y(0)lniif 



This expression shows that the gap vanishes as 
exp(— Ct./(G'(0) — G^)) when the spin-phonon couphng 
constant goes to the critical value, indicating that the 
phase transition between the gapped phase and the 
gapless phase in the antiadiabatic regime is in the 
Berezinskii-Kosterlitz-Thouless (BKT) universality class. 
For fixed G(0), the Eq. H65|l also indicates that there ex- 
ists wbkt, such that for > ujbkt the gap vanishes via 
a BKT transition. The implicit equation giving lubkt 



reads: 



uGjO) 
2a 



lubkt- 



mi 



1 



|y(0)|ln( 



(66) 



which shows that ivbkt is an increasing function of the 
spin-phonon coupling constant. 

The functional dependence of the gap on the spin- 
phonon coupling constant (|65|l is similar to the one ob- 
tained in the case of a frustrated spin chain. This re- 
sult is roughly in agreement with the results of canonical 
transformations that eliminate the phonons from the spin 
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HamiltonianiS*^. Upon closer inspection however, one 
finds that the factor of does not appear in the for- 
mulas giving the spin gap in that case. The reason is that 
the canonical transformations of Refs. 'l3''3l' are valid in 
the limit cjq 3> J and only the instantaneous interactions 
are present, whereas in our theory one needs to renormal- 
ize until the scale ujq reaches J before the interactions can 
be considered instantaneous . Thus, we find that there is 
an intermediate regime, A < < J in which the gap is 
still larger than in the frustrated chain limit. The frus- 
trated chain results are recovered only when loq J. In 
the case of the Heisenberg spin chain, the value of j/(0) 
has been estimated to be of the order of —0.25^ which 
is rather large. For this reason, we cannot make quanti- 
tative predictions on the value of the gap. However, we 
can expect from this large value of the marginal interac- 
tion that the critical value of the spin-phonon interaction 
needed to obtain the spin gap in the antiadiabatic regime 
will be rather large. The behavior of the gap as a func- 
tion of the strength of the coupling constant for the case 
y(0) = as given by Eq. HE20|I is represented on figure |21 
The behavior of the gap as a function of frequency is rep- 



; I — 

0.9 - 

0.8 - 

0.7 - 

0.6 - 

0.5 - 

0.4 - 

03 - 

0.2 - 

01 - 

— 
0.1 



1 10 
uG(0) 

FIG. 3: The behavior of the gap A as a function of phonon 
frequency ujo for a fixed value of G(0) and j/(0) = 0. In the 
adiabatic regime ujo <^ the gap is independent of wo 

and equal to Ao = In the antiadiabatic regime , the 

gap is a rapidly decreasing function of phonon frequency given 
by Eq. 



< I 3° 




FIG. 2: The behavior of the gap A as a function of the spin 
Peierls coupling constant G in the case t/(0) = 0. When 
G S> aujo/u, the system is in the adiabatic regime, and the 
gap A varies linearly with G. When G <C auJo/u, the gap de- 
creases very rapidly with an essential singularity for G(0) — 
described by Eq. 1)65^ . 



whereas, in the antiadiabatic regime, we have: 

,1/2 



(cos2(; 



\ u 



exp 



\ym 



(68) 



An ansatz can be made to describe the crossover between 
the adiabatic and the antiadiabatic regime. In the adia- 
batic regime, G{1) ^ 1 for / such that uiQe'- <^ u/a and 
y{l) <C 1. We can analyze the crossover from the adia- 
batic to the antiadiabatic regime by matching the expres- 
sions obtained in the two cases. This matching procedure 
(see Appendix IeI) predicts that the crossover of the two 
regimes is obtained when uG{G)/a = 2ujo — 2ujc (for the 
case of y(0) = 0). The resulting phase diagram is shown 
in Fig^ where three regimes are visible. As loq increases, 
we go from the gapped adiabatic regime to the gapped 
antiadiabatic regime and finally to the gapless regime. 
As is shown the strength of spin-phonon interaction in- 
creases the size of the gapped regime. 



DISCUSSION 



resented on Fig. |21 Concerning the expectation value of 



S„ • S„+i , we have that 
regime, we find 



(cos2fi 



'cos 2(5 



\ 2TrupLUQ 



-Klo 
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In the adiabatic 



(67) 



We shall here discuss our results for the spin-Peierls 
model and compare with the ones present in the liter- 
ature. The main purpose is to give an overview of the 
informations that can be extracted using the SCHA and 
the RG analysis. We stress that our results are also ap- 
plicable to strongly interacting fermionic systems via the 
Jordan- Wigner transformation. Moreover, for fermionic 
systems, both the case of on-site (Holstein) phonons and 
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FIG. 4; The phase diagram in the {G,ujo) plane, where G is 
the spin-Peierls coupling constant uo the phonon frequency. 
The solid line corresponds to the BKT transition between the 
gapped and the gapless phase described by Eq. 1641 . The 
dash-dotted line is a crossover line between the antiadiabatic 
regime, uG(0)/a ^ ujq and the adiabatic regime uG{0) S> ijJq. 



on-bond phonons can be dealt with as they correspond 
respectively to a term cos2</) or sm2(j) in the action. 



A. Comparison with numerical calculations 

The first result, obtained within SCHA, is the criterion 
for the crossover from the adiabatic to the antiadiabatic 
regime, which is loq ~ Ag, where As is the static gap, and 
not uq ~ J, as could have been naively been expected. 
This criterion was obtained previousl y by a two-cutoff 
renormalization group analysis in Ref. |60| in the case of 
interacting fermions. Such result is in agreement with 
a DMRG study of the XY spin chain coupled to dis- 
persionless phonons of frequency too , where the crossover 
is observed when As ~ i^o"^ Another DMRG study in 
Ref. i21_ also showed that in the case of spinless fermions 
interacting with Holstein phonons, the phase transition 
between the gapful Peierls state (adiabatic) to the gapless 
Luttinger liquid state (antiadiabatic) was also obtained 
when too As. 

The second result, obtained using the RG for the SU(2) 
invariant case concerns the behavior of the spin gap in the 
adiabatic and antiadiabatic cases. We have found that 
for low frequency, Wq ^ As, one is in the adiabatic regime 
with a spin gap given by the mean-field approximation^, 

2 

A ~ As 7^^' ^'-'^ '^o ^ As, in the antiadiabatic 
regime, the spin gap starts to decrease rapidly with the 
frequency. This result is in agreement with the DMRG 
study of Ref. Issl In fact, in this regime the RG analysis 
yields: 



ujQe ^ " exp 



mG(0) 



y{0) 



l-y(0)ln- 



(69) 



fo-^ 2^ > |y(0)|/(l + |2/(0)|In(ziei-V('^o«))), and 
otherwise the gap vanishes. In Ht)9|l . we have assumed 
that y(0) < i.e. that we are dealing with an un- 
frustrated spin chain which would not dimerize sponta- 
neously if g = 0. The three different regimes are illus- 
trated on Fig. 0] An immediate consequence of H69() is 
that for a sufficiently high frequency loq, or for sufficiently 
weak spin phonon coupling constant a BKT transition to 
a gapless state is obtained. This transition is analogous 
to the one that takes place in the frustrated Heisenberg 
chain when J2 < 0.24:Ji^MlMM. 

Our findings for the behavior of the gap, both in the 
adiabatic and in the antiadiabatic regime, are in quali- 
tative agreement with the Quantum Monte Carlo results 
of Ref. [23 on the one-dimensional S — 1/2 Heisenberg 
model. There, it was obtained that for heavy phonon 
(i.e. low-frequency), a static gap was present, while in 
the case of a light phonon (i.e. high frequency), no spin 
gap was observed at the lowest temperature accessible 
in the simulation. It is not obvious whether the ab- 
sence of a dimcrization gap was because the temperature 
was still above the zero temperature gap or because of 
the true absence of a gap above the ground state. In 
either case, these results indicate that the spin gap is 
very strongly decreased with respect to the static result 
when the phonon frequency is increased. Refs. I2^l70l us- 
ing Stochastic Series Expansion methods also found that 
for small spin-phonon coupling and uJo/ J = 1/4 no spin 
gap was obtained, but that increasing the spin phonon 
coupling above a critical Oc = 0.225J caused a phase 
transition from the uniform gapless phase to the dimer- 
ized gap ped phase, in agreement with our results. In 
Ref. I23, the existence of dimerization and spin gap was 
analyzed by Quantum Monte Carlo simulations of the 
Heisenberg model for various spin-phonon couplings and 
phonon frequencies. A phase diagram (Fig. 9 of Ref. l2^ 
was plotted. In agreement with 169(1 . it was shown that 
the critical spin-phonon coupling to induce dimerization 
was an increasing function of the phonon frequency. ED 
studiesiLi^ also show unambiguously the existence of a 
threshold in the spin-phonon interaction to induce a spin 

gap when the phonon frequency is non-zero. A constant 

1/2 

spin gap is obtained when g behaves roughly as ujq 
which seems in agreement with the predictions of 169|l . 
Finally, a similar qualitative agreement is found in the 
one-dimensional Holstein model for spin-^ electrons by 
DMRG study^^ . The system undergoes a quantum phase 
transition between the metallic phase and the Peierls in- 
sulating phase at a finite critical value of the electron- 
phonon interaction. Concerning the BKT universality 
class of the transition predicted by H69|l . in the DMRG 
study of Ref. 23, it was indeed found that as a function 
of the coupling constant the quantum phase transition 
from the gapless to the gapped state is a KT transition. 
A KT transition was also found in the related ferniionic 
case in Refs. l2lll65l Another result of Ref. I23I regards the 
evaluation, by finite-size scaling, of Luttinger exponent 
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Kp of the spinless fermions. It was found that approach- 
ing the transition at a critical value of g^^ Kp has small 
deviation from 1/2 from 0.59 to 0.42. Such finding is in 
agreement with a BKT transition driven by the operator 
cos 40 since the value of the Luttinger exponent at the 
transition is then K = 1/2. 

In the antiadiabatic regime we predict a power-law re- 
lationship between the critical spin-phonon coupling gc 
and the frequency wq- We note that a power law relation 
between the critical spin phonon coupling was found in 
Ref. I23. However, the results of Ref. |23 are obtained in 
the limit of XY anisotropy, so that a direct comparison 
of the exponents is not possible. We can make a more 
direct comparison with the data of Ref. If we call 
go A the spin phonon coupling constant used in Ref. llSL 
it is related to our spin-phonon coupling constant by: 



9 



(70) 



yielding for the dimensionless spin-Peierls coupling con- 
stant: 



G 



'Jda 



Tt'^LOq 



(71) 



The exact diagonalizations in Ref. [T^ were performed for 
ujQ — 0.3 J. The dimensionless parameter u/{aujQ) — 
^ ~ 5.2, indicating that a continuum description such 
as ours should be still applicable. The values given in 
Ref.0, Fig. 1, lead to a dimensionless coupling constant 
in the range [0.17, 0.27], which is at the limit of the per- 
turbative regime. Since TrG/2 < .429 < 2luq/J = 0.6, 
we are in the antiadiabatic regime, not far from the 
crossover. The gap we are calculating is since our 
RG approach does not take solitons into account. The 
result of the comparison is shown on Fig. |31 where we 
have replotted the data of Ref. fii for the gap A'^^ , along 
with the formula (jE20(l in which logarithmic corrections 
are neglected. Obviously, the overall behaviors of the 
gap with the coupling constant are very similar. On a 
more quantitative level there is a discrepancy between 
our results and those of Ref. by roughly a factor 2. 
Such difference is clearly not due to the logarithmic cor- 
rections. Because our coupling constant is already rather 
large, the logarithmic corrections should be rather small. 
Moreover, if the discrepancy was caused by logarithmic 
corrections, it would diminish as the coupling constant 
increases, which is not the case here. This leaves us with 
three possible explanations of the quantitative difference 
between our results and those of Ref. Ts'. The most likely 
explanation is that a coupling constant G in the range 
[0.17, 0.27] is already a rather large value of the coupling 
constant, and a one-loop RG such as ours is not sufficient 
to obtain the gap quantitatively in this regime. A two 
loops or higher order calculation may reduce the scale 
of the gap with respect to one loop and lead to better 
agreement. However, since the usual techniques for deriv- 
ing renormalization group equations^i^iSiS are restricted 
to one-loop, one needs to develop a field theoretic RG 
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FIG. 5: Comparison of our prediction for A''^ with the results 
of exact diagonalization(Ref. IT^ . In both cases, the gap is 
an increasing function of the spin-phonon coupling constant. 
Our results are roughly twice the exact diagonalization result. 



approacb24 to check this. Such an approach is beyond 
the scope of the present paper. Alternatively numerical 
data for a smaller coupling constant would be interesting 
to compare to. The second possible explanation is that 
our procedure to match the results of the RG in the adi- 
abatic and the antiadiabatic is introducing an incorrect 
scale factor in the antiadiabatic limit. This is possible 
if for instance, when we are in still the adiabatic regime 
but near the crossover, the interaction with the phonons 
is causing a reduction of the gap. Such an effect is ignored 
in the Eq. HF20|I . A possible last explanation is that be- 
cause the ratio of the phonon frequency to the exchange 
coupling is still not very small, the continuum treatment 
is not sufhciently accurate. Note however that despite 
the relatively extreme case of this numerical data (large 
coupling constant and phonon frequency) with respect to 
a continuum approach and first order RG analysis, the 
quantitative agreement is still reasonably good. 

In addition of providing an analytical framework to 
describe the behavior of the spin Peierls gap as a func- 
tion of coupling constant or frequency our analysis al- 
lows to extract other physical quantities. In particular 
a quantity that can be deduced from our calculations 
to describe the Peierls ordering structure of the ground 
state is the dimerization S = (qn)- We can calculate it 
from the magnetic order parameter through the relation. 



(-)"-f (cos(2(/))). In the adia- 



batic regime, the results of the RG analysis gives: 



1/3 



whereas in the antiadiabatic regime we obtain: 



(72) 



exp 
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(73) 
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Since by scaling, the expectation value of (cos(20)) is re- 
lated to the spin-Peierls gap by the relation (cos2(/)) ~ 
(y)^, we immediately obtain the scaling relation for the 
dimerization order parameter, 5 ~ ^(■7)^- It was found 
in ED studieSiiS [Fig.l of Ref.[ll| that the alternation S 
in the exchange integral and the spin gap were increas- 
ing functions of the spin phonon coupling constant with a 
threshold. The correlation of displacement and dimeriza- 
tion predicted by J7^ and (|75|l are also neatly illustrated 
on Fig. 8 of Ref. IH"^ 

From the behavior of (cos(20)) in the adiabatic and 
antiadiabatic regime we can also immediately infer some 
features of the spin-spin correlations function. In the 
antiadiabatic regime, from H68|l . we deduce that for a 
fixed frequency, the correlations decrease at increasing 
the spin- phonon coupling [Fig. 15 of Refl2^. At a fixed 
spin- phonon coupling and increasing uo, the spin-spin 
correlation functions become instead less and less affected 
by g. These findings are again in agreement with the re- 
sults obtained in Ref.|2^ We thus see that our derivation 
provides a unified framework explaining and generalizing 
the previous studies. 



B. Relation to experiments 

Let us now turn to experimental systems. Our re- 
sults clearly show that non-adiabatic phonon dynamics 
strongly renormalizes the magnetic correlations and the 
dependence of the gap. 

A well known example of a material where such 
strong renormalizations are observed is the spin-Peierls 
materials CuGeOa. In this material, the phonon fre- 
quency is rather high compared to the actual spin gap 
{ujQ ~ 310K). Interestingly, the thermodynamics of 
this material can be fitted with a frustrated spin chain 
model2^, with J2/J1 = 0.36 i.e. well into the spin gap 
regime. As was pointed outi^ii^ such a dimerization is 
not intrinsic but due to the spin-phonon coupling itself. 
Indeed our RG analysis shows that the low energy prop- 
erties of a spin chain coupled to dynamical phonons are 
similar at low energy to those of a frustrated spin chain 
provided that the phonon frequency is above the zero 
temperature spin gap, which is the case in CuGeOa. In 
such a system one can expect a very strong reduction of 
the gap due to the finite phonon frequency as shown in 
Fig. 13 and was noted before for CuGeOa. As a point 
of comparison of the interplay between the phonon fre- 
quency and the phonon coupling constant we have re- 
ported on Fig.l^the various compounds listed in Table 1 
of Ref. Ha. From Eq. gij), we expect Amf ^ ^G(O), 

and therefore using Eq. (|E20p . A = ujQe^^'^^ . Note 
that the values of the gap taken here are only indica- 
tive, in connection with our one-dimensional analysis. In- 
deed, they are (i) dependent on the measurement method 
and slightly differ depending which quantity is measured, 
(ii) are dependent in part of the interchain couplings 
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FIG. 6: Various compounds exhibiting a spin-Peierls transi- 
tion plotted on a gap A/wo vs Cross-Fisher gap Amf/u>o- 
The data for the gap A, the phonon frequency ujq and the 
Cross- Fisher gap Amf are extracted from Table 1 of Ref. 17^ 
Note that the experimental determination of the gap is to be 
taken with a grain of salt given the differences between the 
various determinations from various measurements (neutrons, 
thermodynamics, etc.) The full line is Eq. IIE20II which is the 
dependence of the gap on frequency that we expect when log- 
arithmic corrections to scaling can be neglected. 

which we have not treated in the present theory. We 
nevertheless see on Fig. that the agreement between 
our calculated values of the spin-Peierls gap and the 
observed one for various systems both in the antiadi- 
abatic regimes is quite decent. Although CuGeOa is 
the material for which the effects of the phonon fre- 
quency are the stronger, another material for which the 
present study could be relevant is MEM(TCNQ)2. For 
this material the phonon frequency is of the same or- 
der than the spin-Peierls gap A '--^ 30 — 60K, and one 
can thus still expects effects of the finite phonon fre- 
quency on the spin-Peierls gap. The other two com- 
pounds (TMTTF)2PF6 and (BCPTTF)2PF6 are closer 
to the adiabatic regime (since they have A^^^ > ojq) 
and thus are expected to have a gap less dependent of 
the phonon frequency than CuGeOa or MEM(TCNQ)2. 
We note that for these two compounds, the agreement 
with our formula is not good. For (TMTTF)2PF6, this 
might result from the fact that the charge localization 
temperature Tp ~ 200 K is relatively low compared to 
the spin-Peierls transition temperature and charge fluc- 
tuations can still influence the transition. Indeed, a de- 
scription based on adiabatic phonons interacting with 
both charge and spin fluctuations can successfully ac- 
count for both the magnetic susceptibility and the NMR 
relaxation rate in this material^. A related explanation 
of the discrepancy could be the existence of a charge or- 
dering transition^6,77 (TMTTF)2PF6. Such a tran- 
sition can affect the mean-field spin-Peierls transition 
temperature of the material2^*22i and thus invalidate our 
simple minded estimate of the spin-Peierls coupling con- 
stant. Finally, in both TTF materials, antiferromagnetic 
interchain coupling could be relatively important, and 
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they may diminish the spin-Peierls ordering resulting in 
a smaller spin-Peierls gap**^. 

In order to test for further prediction of the above 
determination of the gap it would be interesting to be 
able to vary continuously the phonon frequency. Pres- 
sure could be an interesting way to address this question. 
Since when applying pressure both the exchange constant 
and the phonon frequency are to be affected, one has to 
compute the net effect on the gap which our theory al- 
lows to do. Such measurements could allow to follow the 
behavior of the gap such as described in Fig. |3| 



VI. CONCLUSIONS 

In the present paper, we have analyzed spin-Peierls 
problem for a single spin- 1/2 chain coupled to an op- 
tical phonon of frequency loq using bosonization tech- 
niques. The bosonized action was approximately solved 
by using the self-consistent harmonic approximation^iiS 
In the low frequency limit, we have reproduced the re- 
sult obtained by Cross and Fisher by the mean-field 
approximation^ In the high frequency limit, we have 
shown that the retardated interaction was giving rise 
to a term that was local in time and was identical to 
the bosonized form of a frustrating next-nearest neigh- 
bor interaction in agreement with the canonical trans- 
formation approach. The self consistent approxima- 
tion also allowed us to describe entirely the crossover 
between the two regimes, regime that we have con- 
sidered, the crossover frequency was given by the spin- 
Peierls gap As calculated for a static phonon (w — > 
with fee — niLUQ fixed). The adiabatic regime extends in 
the region luq <^ As, and the antiadiabatic regime ex- 
tends in the region ujq ^ Ag. All the previous findings 
can be recovered by the renormahzation group by us- 
ing a two-step approximation as in Ref. [SO- We stress 
that although we use the same two-step approximation 
as Ref. 01 the behavior of the spin gap that we obtain 
in the low frequency limit in our renormahzation group 
is different from the one that one would deduce from the 
renormahzation group applied to spinful fermiqns at half- 
filling coupled to optical phonons as in Ref. 60. The rea- 
son for this is that in our problem the charge mode is ab- 
sent, making the spin-phonon interaction more relevant 
than in Ref. l6Ct The advantage of the renormahzation 
group approach over the selfconsistent approximation is 
that the former is applicable in the SU(2) invariant case 
where the induced non-retarded approximation becomes 
marginal. For the SU(2) invariant case, the crossover 
frequency between the adiabatic and the antiadiabatic 
limit remains wq ~ A^. However, due to the marginality 
of the induced term, a BKT transition in the antiadi- 
abatic limit between the gapped state and the gapless 
state becomes possible. Near the transition, the spin 
gap A drops very rapidly with the phonon frequency as 
A e "^c-^o and vanishes for luq > Uc- These results 
are in qualitative agreement with numerical studies^ and 



with the canonical transformation method. ^-'^ As the fre- 
quency LJc is a increasing function of the spin-phonon 
coupling, for fixed loq there exists a critical spin-phonon 
coupling below which the spin gap disappears. We have 
also examined in the light of the present theory the ex- 
isting experimental compounds exhibiting a spin-Peierls 
transition as summarized for example in Ref.ll2l We find 
a good qualitative agreement with the dependence of the 
gap predicted by our theory. 

Our analysis thus provides a unified analytical frame- 
work in which to analyze the spin-Peierls transition. It 
leads the way to interesting extensions. In particular on 
could expect to tackle with similar methods the effects 
of the interchain couplings, or the effect of impurities on 
the spin-Peierls transition (see e.g Ref. [jl! and references 
therein) . 
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APPENDIX A: EXPRESSION OF THE 
PROPAGATOR IN THE SCHA 

We need to obtain an expression for the propagator G. 
We have: 



G{x,t) 



t:Ki 
ttK 



E 



dq 



i{qx- 



r) 



27r Lul + u^{q^ +^~^) 



G{x,t) 



(Al) 



The reduced propagator Q satisfies to the partial differ- 
ential equation: 



g{x,T) = -S{x) J2 SiT-nP){A2) 



n— — oo 



and to the following properties: 

g{x,T + p) = g{x,T) 
g{±x,±T) = gix,T). 



(A3) 



To solve the partial differential equation ljA2|l . we con- 
sider first the auxiliary partial differential equation: 



goix,T) ^ -S{x)d{T). 



(A4) 



16 



Clearly, if we have a solution of (jA4|) , we can easily de- 
duce from it a solution of l|A2|l : 



G{x,t) = ^ go{x,T-n(3). 



(A5) 



n— — oo 



An explicit solution of the PDE (|A4p is readily found by 
Fourier transformation, and application of Eq. (9.6.21) 
of Ref.EE One has: 



goM = -Koi — - 



(A6) 



and thus: 



We note that in reality we have slightly cheated. Be- 
cause of the cutoff on the momentum integral in l|Aip . 
the function S{x) is in fact smeared out into a function 
d\{x) which goes to the delta function only for A —^ oo. 
Consequently, the true propagator Qx is in fact the con- 
volution of S\{x) with the function defined by Eq. IjATp . 
A simple way to incorporate the cutoff in (|A7(I is to per- 



form the replacement x 



with \-^. This 



substitution has the advantage of preserving the symme- 
tries Finally, we have: 



^ x^ + u2(t — -|- o? 



(A8) 



27r ^ " \ i 

n—~-oo \ 



It is easily seen that the series in (|A7|) is convergent, and 
that the function defined by (|A7|) satisfies all the con- 
ditions \K?t\ . The result (|A7(I could also have been ob- 
tained by using the Fourier transform of the Dirac comb. 



APPENDIX B: SOLUTION OF THE 
VARIATIONAL EQUATIONS 

In the present section, we study the solution of the 
variational equations derived from minimization of the 
free energy H41(l with respect to ^. The contribution of 
the region r ^ ^/u can be written as: 



4(7ra)2p(jj2 

^2 



e 'a 
"2^ 



2K 



2e '^ujq^/u 



9 





2 A' 



2K 



I e'u 



2K 



7(1 



+ 2X,2e-^^) 



-2K 



'2K 



-2K 



where 7 is the incomplete gamma function (see Ref. ISS 
chap. 6 p. 260). Using the identity: 



The contribution of the region t 3> ^ can be rewritten 
as: 



, 2K 

a 



u 

2K 



-2K 



) 7(l-2ii:,2e" 



(B2) 



r(l - 2K) - r(l - 2K, 2e-'^ — ) 

u 



and noting that the first term in the right hand side is 
independent of ^, we can rewrite up to a renormalization 
the short-distance contribution to the variational free en- 
ergy as: 



, 9 I 

'Ai^afpul \ 2C 

9' 



2K 



dve-" 



2{-KafmLjl V 2^ 



2K 



(B4) 



2K 



7(l-f2if, 2e-^^) 
u 



2^ 



r(l - 2K, 2e- 



(B3) Thus, the variational free energy to use reads: 
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F = Fo- 



V u 



2K 



4(7ra)2p(jjQ 
T{1-2K, 2e^^— ; 



2e 



2K 



-2e" 



(B5) 



APPENDIX C: DERIVATION OF THE 
RENORMALIZATION GROUP EQUATIONS 

To derive the renormalization group equations, we 
start from the action H52I) where the function ^^^^^{t) 
is defined in ((T^ and satisfies: 



Renormahzation group equations are obtained from Op- 
erator Product Expansion techniques (OPE)j22iS The 
following OPEs are needed: 



/3/2 



13/2 



T^uo, l3{T)dT = 1 



(CI) 



cos 2(j){x, t) cos 2<j){x' , t') 



cos 40(2;, r) cos40(x', r') 



1 - i (2(x - x')d,(l){x, t) + 2(r - T')drc^{x, T)f + cos4(/.(x, r) 



cos2(/)(x, r) cos4(/)(x', r') ~ - cos 2(/)(x, r) 



1 - ^ (4(x - x')d^c^{x, r) + 4(r - t')9.0(x, t))' 



(C2) 
(C3) 
(C4) 



To find the renormalization group equations, we write the 
partition function, and expand to second order around 
the Gaussian fixed point. The important contributions 
are of first order in (due to the nonlocality of the action 
(|52|l ). and of second order in g^^ and g±g'^. Then, we 
change the cutoff by a ^ ae'^' . We then apply the OPEs 



(|C2|) to obtain the short distance contributions of the 
terms with < {x — x')^ + w^(r — r')^ < a^e^**' to the 
renormalization of the coupling constants. Proceeding in 
that way, we obtain the following corrections 0{dl) to 
the action: 



2/9W, 

coming from the g^ term. 



(—Y T^^oWi-) I dxdrdr'i 



-(r-T')n9r<^r + 7rcos40]. 



(C5) 



1 f ^9A 



,1 



2u \(27rCk)2y JaKrKae" 

coming from the g±^ term, and: 



rdrdO- 



1 - 8r2 cos2 e{d^(l)f - 8— sin^ e(drC, 



(C6) 



-25± „ 1/5 



(27ra)2 2pw, 



dxdrdr' cos2(t){x,T')V^^m{T - t') X 2 dx" dr" - cos 2<j){x,T) (C7) 



a<|r — r"|<ae^ 



coming from the g±g'^ term. The term (jCSII contributes tributes to the renormalization of u, K. Finally, the term 
to the renormalization of K,u,g±. The term (|C6|) con- 
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r,dl 



(jC7|) contributes to the renormalization of . After hav- 
ing performed the mode integration, we have to restore 
the original cutofS^. An operator O of scahng dimen- 
sion do is rescaled as O ^ Oe^'*°' whereas coordinates 
are rescaled as : 
(3 /Je-'" i.e. 
tion, we obtain g± — > gj^e^^"*^)'*'. The rescaling of 
is a bit more subtle. First, we notice that the rescaling 
of a — > ae'^'' amounts to a rescaling of ujo,l3 inside the 
function V. Hence, the rescaling acts on V as: 



xe"" and r re . One also finds 
Te'^K As a result of this opera- 



V, 



V. 



dl 



(C8) 



APPENDIX D: LOGARITHMIC CORRECTIONS 
TO SCALING 

Let us consider a spin-1/2 chain with a static dimer- 
ization, described by the Hamiltonian: 



i/ = J2^(l + <5(-)")S„-S„+i, 



(Dl) 



Bosonization and scaling arguments2i4i lead to the pre- 
diction of a gap A ~ S^^^. However, the presence of 
a marginally irrelevant operator induces corrections to 
scaling^2ii2iSiSiS4 and the gap behaves in fact as: 



and this way of rescaling T^^jo guarantees that the con- 
straint (jCl|) remains satisfied. The rescaling of V ab- 
sorbs the opposite rescaling of one of the components 
x,t,t'. As a result, the rescahng of g^ is given by 
g"^ g(2-2A:)dz^2^ notice that this result is in con- 
trast with the case of a disordered system^ in which the 
disorder D is rescaled as I? ^ £)^(3-2K)di _ Mathemati- 
cally, the diff'erence arises because in the case of phonons, 
we want to keep the weight in the function V constant 
under the RG flow. This constraint means that physi- 
cally we are converting the non-local interaction that ex- 
ists at high energy into a local one represented by the g± 
term by integrating out successively its short distance 
contributions. The remaining weight is then given by 



2 f 



/3/2e 



V, 



{T)dT and goes to zero as I is increased. 



Adding all the 0{dl) contributions, both from integration 
of short distance terms and from rescaling, we finally ob- 
tain the renormalization group equations as: 



— ( ^ ^ ^ ( 3 V 1-1 



at \27rwA / u \iTu/ pujQ Zn^u^ 
d_ f u \ ^ 
dl \27rK) 

dl \TTU/ TTU VtTU/ pUjQ 

dl yiTupujQ J 
dujQ 

nr 

dT 

H 

(C9) 



nu/ nupujQ 



T 



We note that for g = these equations reduce to the 
usual RG equations of the sine-Gordon modelS. If we 
neglect the contribution of g± to the renormalization of 
g^, we see that the gap obtained for g"^ /{upoj^) ~ 1 coin- 
cides with the gap predicted by Cross and Fisher^. Tak- 
ing g± into account for K ^ 1/2 leads to logarithmic 
corrections in the dependence of the gap on G. These 
corrections are discussed in Appendix IdI 



A = 



1.723,52/3 



1/2 



(l + ll 



(D2) 



As a result of logarithmic corrections, the ground state 
energy of the dimcrizcd spin chain behaves as: 



T 



0.2728^4/3 



y2(0)|ln 



ly(o) 



1.36125 



(D3) 



If we now consider that the chain ljDl|l is coupled to 
adiabatic phonons, one needs to minimize the total en- 
ergy ke/g^5^ — B5^/^\\n5\~'^ with respect to 5 yielding 
(52/3| ln(5| ~ g^ /ke- The gap thus behaves as 



A = 6.07- 



In- 



-3/2 



(D4) 



where we have used the factors quoted in Eqs. (|D2j|- HD3|) 
and Ref4^. We know show how the expression (|D4() can 
be recovered within our RG approach. Using the initial 
conditions with SU (2) symmetry, the RG equations are 
given by the Eqs. H61I62|I 

If we now assume that G(0) ^ 2/(0), in Eq. 1)61(1 we 
can take G = 0, so that the previous equations reduce to 
the single BKT equation 



dy_ ^ 2 
dl ^ 



(D5) 



For y < Q, this equation fiows to a fixed point y* = 
with the flow given by: 



yii) 



y(o) 



l-y{0)l 



(D6) 



Using HD6() in Eq. (|62|) we can easily integrate it and 
obtain: 



G(0 = G(0) 



[i + |y(o)|;]3/2- 



(D7) 



This equation should break down for G(/o) — |y(^o)l £^ud 
after that an exponential flow of G is expected. Using 
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Eq. ljD7|l , the strong coupling behavior is obtained when: 



G(0) 



(D8) 



(D9) 



Solving Eq. l|D9p by iteration and using the first iteration, 
one finds the following scaling for the spin-Peierls gap 
valid for small G: 



G(/o)e-'» 
where Iq is given by: 

G(Oy« = |y(0)|(l + |z/(0)M^/'. 



A, 



G(0) 



1 



|y(0)|lnM^ 



3/2 ■ 



(DIO) 



This behavior is in agreement with the behavior obtained 
in Eq. ljD4|) by considering logarithmic corrections to the 
energy of the dimerized chain. 



APPENDIX E: STUDY OF THE SU(2) 
INVARIANT RG EQUATIONS 

We consider the system of 2 coupled first order dif- 
ferential equations (|61|l - l|()2|l . A convenient approach in 
the antiadiabatic regime is to recast this differential sys- 
tem as a single differential equation, ft is also convenient 
to make the approximation G{1) ~ G(0)e' as to render 
second order equation linear. Introducing^: 



Y{1) = exp 



yil')dl' 



(El) 



We obtain the following second order differential equa- 
tion: 



2u 

uG{0) f uJoae 
2ujQa 



(E2) 



Yil) (E3) 



with the initial conditions for 1^(0) = 1 and Y'{0) = 
-2/(0). 

It is possible to simplify the second order differential 
equation (|E2|) by a variable change to s = Writ- 
ing Y{1) = Z{^e^), one finds that Z{s) satisfies to the 
differential equation: 



d^Z IdZ uG{0) 



ds^ s ds 
with initial conditions: 



LOoa ^ 



e-^Zis) = 0, 



LOoa dZ ,ujQa ^ 
u ds u ' 



-y(o). 



(E4) 

(E5) 
(E6) 



The equation (|E4|I can be recast in the form of an integral 
equation which reads: 



Z{a) = l-y(0)ln 



uG(0) r ds' I-'" 



s"e-' Z{s")ds"{Yn) 



We study first the antiadiabatic limit, MG(0)/a ^ loq. 
Let us begin with the case 2/(0) = 0. By iterating the 
equation IIE7p once, we obtain: 



Z{s) = 1- 



uG(0) 
2^0 



V U I \ LOnOL 



V . ) 



.o(f^), 
2wo 



(E8) 



where Ei is defined in Ref. 58. Using the fact that 
u! — Oa/u <C 1, and ujQa/ue'' 3> 1, we can rewrite this 
equation as: 



Y{1) = 1 



iG(0) 
2wo 



? - In ( — 1 + 7 - 1 

ujQa 



(E9) 



where 7 = 0.577 ... is Euler-Mascheroni's constant. 
The gap is obtained when Y{1*) — 0, i.e. 



H 7^r7T-r + 1 - 7 



\aLUo J uG{0) 



(ElO) 



The differential equation (|E4p can also give some in- 
dication on the crossover scale to the adiabatic regime. 
Assuming that -^^^ ^ 1, it is reasonable to replace the 



term e * with 1 in (|E4p , yielding the approximate equa- 
tion: 

ds s ds ZujQa 

The solution of the above differential equation is easily 
found^^ in terms of Bessel functions. One has: 



(E12) 

The initial conditions yield the following linear system 
for A, B: 



AJn 



AJi 



woaG(0)V/^ 



2u 

2ir 



BYn 



BYi 



' f ojQaG{Q )\^'^'' 
V 2^ 



2u 



woaG(O) 



= 1 

- y(o) 

(E13) 
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If we consider the case of y{0) = 0, we find: 

For ujQaG{0)/u ^ 1, an approximate solution is Z(s) — 
Jo((uG(0)/(2woa))i/2s). We obtain Z{s) = for 
{uG{0)/{2ujoa)y/^s = jo,i - 2.40482 (cf. Ref-IH). The 
resulting strong coupling scale would be given by: 





(-)l 


[g(o) 


\uioaJ _ 



(E16) 



We note that at this scale G{r*) > 1, but ^e'* < 1. 
This indicates that G{1) is reaching strong coupling be- 
fore the scale I* is reached, and in this regime, the gap is 
produced by G directly. The true strong coupling scale in 
this adiabatic regime is then l*^ = ln(l/G(0)). The above 
calculation thus allows to determine when the crossover 
from adiabatic to antiadiabatic regime is obtained. 

Comparing Z*^^ with l*ntiad ^^'^ ^^^^^ their difference 
is minimal when -^^^j^ = 2. Thus, for mG(0) <^ ujoa, we 
are in the antiadiabatic regime, and: 

A^LUoe^^^'^'^e'^^ (E17) 

For wG(0) > LUoa, the gap behaves as A = CuG(O). We 
need to match the two results. To do this, we require that 



the two expression of the gap are equal and have the same 
derivative with respect to G(0) when uG(0)/q! = 2ujq. 
This yields the following ansatz for the gap: 



A=- — G(0),-G(0) > Wo (E18) 
2e a a 

A = Cujoe , -G(0) < loq (E19) 
a 

The constant C can be obtained using results for the adi- 
abatic limit. Using Eq. (8) of Ref. [s^ we find that 
C/(27re) = 0.627 leading to C = 10.7. This leads us to 
the following ansatz for the gap: 



A = 0.627^^ for 2wo < -G(0) (E20) 

A = 10.7woe"^''^ for 2wo > -G(0) (E21) 

a 

The above ansatz Eq. HE20|I tends to overestimate the 
the gap in the adiabatic regime as it neglects completely 
the effect of the non-zero frequency on the gap. However, 
the lack of a precise criterion to decide when the RG flow 
has reached the strong coupling regime prevents us from 
finding a better answer. 

We now turn to the case of y(0) ^ 0. We first look 
at the antiadiabatic limit. We obtain by iterating the 
equation HE7|I that: 



Z{s) = l-y(0)ln( — 
\uJoa 



uGjO) 
2uJoa 



1 + ^) e-^ - ,(0) (e-^ + E, i^))] In ( ^ 
u / V \ u //J \uJoay 



mG(0) 
2u}oa 
uG{0) 
2ujaa 



1 _ j^(0) In 



2/(0) 



2E, 



uJoa 



2G 



In 



E, (—] In ^Ei (s) In s - 2Ei (s) - 2G(s) - e"^ In s 
\ u / u 



(E22) 



Using the expansions of Ei and G for small argument: We find the following expression for Y(l) at small 

uG(0)/(woa): 



E^{s) = --i-\ns + o{l) (E23) 

(E24) 
(E25) 



G(s) = -i(lnsf + ^^ + 0(1) (E24) 
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Yil) 



2uJoa 

-r"(i 



1 - 7 + In 



27] 



u 
Wo a 
mG(0) 



£2/(0) 



2tJoCk 



In" 



+ 2(1 -7) In 



7 + ln( — ) 



r 



For y(0) < a solution exists if: 



uGjO) 
2uJoa 



> 



m\ 



(E26) 



l2/(0)|ln(^) 
In other words, the existence of the gap is controlled by 



the ratio of the coupling constant to the marginally irrel- 
evant perturbation measured at the energy scale uiQe^^'' . 
We haveS r"(l) =7-^ + ^. However, for what follows, it 
is convenient to make an approximation r"(l) ^ 7^ + 1. 
Then, one can rewrite: 



Y{1) = 1 



uGjO) 



1 - y(0) In 



=1-7 



In 



V 



y(o) 



uG(0) 



l-j/(0)ln 



ue 



1-7 



(E27) 



The strong coupling scale is obtained for Y{1*) = which 
gives us: 



I* = 1-7 + ln 



1 



uG{0) 



y{o) 



(E28) 



2woQ i_y(o)ln(iifiljl) 

in the above formula gives 



We note that making y(0) 
back the expression IjElOjl . 

Turning to the crossover to the adiabatic regime, we 
note that again, due to the lack of a precise criterion for 
cutting the RG flow, we can only propose an ansatz to 
relate the two regimes. Extending the reasoning made in 
the previous discussion of the case y{0) = 0, we expect 
that the crossover happens when l*^ — l*nuad minimal. 
The lengthscale l^^ has been obtained in section IdI It 
reads: 



= - In G(0) + - ln[l - |y (0) I In G(0)] (E29) 
Minimizing the difference of lengthscales then gives: 

3 \ym 



uG{0) 
2uJoa 



1 



\ym 



uGjO) _ 
2uJoa l+|y(0)|ln- 



21-|2/(0)|lnG(0) 



(E30) 

In the case of ?;(0) <C 1, one can write ^G(O) — 2u!o{l+£). 
The quantity e is straightforward to obtain by expanding 
(|E30|I . Finally, one has: 



-G(0) 



2wo 



mi 



\y[ 



l2/(0)|ln- 



2 1 + |y(0)|ln 



(E31) 



Therefore, the crossover scale is only weakly affected by 
the presence of logarithmic corrections. An ansatz sim- 
ilar to Eq. (jE20|) can be derived in the case of a spin 
chain with marginally irrelevant operator. Using the re- 
sults of Refi^, combined with the Adp.IdI we obtain the 
following expression of the gap in the adiabatic regime. 



A = 1.96886 -G(0) 
a 



2/3 

557G(0) 



(E32) 



Our calculation of the crossover scale shows that this 
expression is valid provided 2uJo < uG{0)/a. For 2a>o > 
uG{Q)/a, an expression of the gap of the form A ~ Ce~' 
with r given by Eq. HE28P is valid. Matching the two 
expressions for 2wo = uG(0)/a yields: 



A = 10.704 t^o- 



X exp 



-1 



uG(0) 



\ym 



2'^"" l+|y(0)|l„(^) 



for 2cjo >< uG{0)/a. Letting y{0) 0, in Eqs. 
(|E33(I . we recover the previous formulas (|E20(I . 



(E33) 



(IE321) 
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